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ABSTRACT 

We model the number of detectable supernovae (SNe) as a function of redshift at different flux thresholds, 
making use of the observed properties of local SNe, such as their lightcurves, fiducial spectra, and peak magnitude 
distributions. We assume that the star-formation rate (SFR) at high redshift traces the formation rate of dark matter 
halos. We obtain a rate of 0.4-2.3 SNe/arcmin 2 /yr at z ^5 5 at the near infrared (4.5 /im) flux density threshold 
of 3 nJy (achievable with the James Webb Space Telescope (JWST) in a 10 5 sec integration). In a hypothetical 
one-year survey, it should be possible to detect up to several thousand SNe per unit redshift at z ~ 6. We discuss 
the possible application of such a large sample of distant SNe as a probe of the epoch of reionization. By heating 
the intergalactic medium (IGM) and raising the cosmological Jeans mass, the process of reionization can suppress 
star formation in low-mass galaxies. This could have produced a relatively sharp drop in the SNR around the 
redshift of reionization (z re ). We quantify the detectability of this feature in future surveys of distant SNe by 
varying the redshift and duration of the feature, as well as its impact on the SFR in low-mass halos, which results 
in different redshifts, widths, and sizes of the drop in the expected SFR. We find that the drop can be detected out 
to Zm ~ 13, as long as (i) the reionization history contains a relatively rapid feature which is synchronized over 
different regions to within Az < 1 - 3, (ii) the star-formation efficiency in halos that dominate reionization is 
e* ~ 10%, and (iii) reionization significantly suppresses the star formation in low-mass halos. Depending on the 
details of (i) - (iii), this could be achieved with a survey lasting less than two weeks. Detecting this signature 
would also help elucidate the feedback mechanism that regulates reionization. 

Subject headings: cosmology: theory - early Universe - supernovae: general - galaxies: high-redshift - 
evolution 



1. introduction 

The formation of the first astrophysical objects and the sub- 
sequent epoch of reionization offer a wealth of insight into the 
physical processes of the early universe. Two recent develop- 
ments, revealing the impact of reionization on the Lyman a ab- 
sorption spectra of z k; 6 quasars (Fan et al. 2004; Mesinger 
& Haiman 2004; Wyithe & Loeb 2004b), and on the tempera- 
ture/polarization cross-correlation anisotropy in the cosmic mi- 
crowave background (CMB) radiation (Kogut et al. 2003; Spergel 
et al. 2003), suggest that reionization history may have been 
very complex (see, e.g., Haiman 2003 for a review). Various 
tools designed to probe this epoch have been proposed, relying 
on upcoming 21-cm surveys (e.g. Madau et al. 1997), high- 
redshift quasar and galaxy spectra (e.g. Mesinger et al. 2004), 
high-redshift galaxy surveys (e.g. Barkana & Loeb 2000), GRBs 
(e.g. Barkana & Loeb 2004a), and more. Nevertheless, there is 
currently a striking dearth of data or theoretical consensus con- 
cerning this important period in the universe's evolution. Due 
to the considerable technical challenges associated with such 
high-redshift observations, some of the most useful objects in 
the near future are probably going to be some of the brightest 
and hence easiest to detect, such as bright quasars, GRBs and 
SNe. Here we make use of the observed properties of local SNe 
to construct detailed high-redshift SNe rates, obtainable from a 
future infrared space telescope such as JWST. We also show that 
such a sample of high-redshift SNe could be an invaluable tool 
in probing the reionization epoch, and the physical feedback 
mechanism that regulates it. 

There have been several previous studies of the expected early 
SNR. However, these either did not focus on the expected rates 
at redshifts beyond reionization (e.g Madau et al. 1998), fo- 



cused exclusively on very high-z rates from the first generation 
of metal-free stars (Wise & Abel 2005), and/or did not address 
the impact of reionization on the SFR (e.g. Miralda-Escude & 
Rees 1997; Dahlen & Fransson 1999). Likewise, previous the- 
oretical predictions of the high-redshift SFR (Barkana & Loeb 
2000; Bromm & Loeb 2002; Choudhury & Srianand 2002) as- 
sumed a fixed degree of suppression due to reionization, match- 
ing numerical simulations which did not include self-shielding 
(Thoul & Weinberg 1996). The main distinctions of the present 
paper are that we compute the expected SNe rates for the wide 
range of redshifts within which reionization is expected to have 
occurred (6 ^ z re ^ 20), and we allow various degrees of photo- 
heating feedback. We also supplement the standard halo-based 
estimates of the SFR with a more elaborate estimate of the cor- 
responding observable SNR, utilizing the observed properties 
of low-z SNe. 

In the presence of an ionizing background radiation, the IGM 
is photo-heated to a temperature of <; 10,000K, raising the 
cosmological Jeans mass, which could suppress gas accretion 
onto small-mass halos (e.g. Efstathiou 1992; Thoul & Wein- 
berg 1996; Gnedin 2000; Shapiro et al. 1994). Reionization is 
then expected to be accompanied by a drop in the global SFR, 
corresponding to a suppression of star formation in small halos 
(i.e. those with virial temperatures below T wu <, 10 4 - 10 s K). 
The size of the drop is uncertain, since the ability of halos to 
self-shield against the ionizing radiation is poorly constrained 
at high redshifts. Early work on this subject (Thoul & Weinberg 
1996) suggested that an ionizing background would completely 
suppress star formation in low-redshift "dwarf galaxy" halos 
with circular velocities v c ; lc <J 35 km s" 1 , and partially suppress 
star-formation in halos with 35 km s -1 <, v c ; rc <, 100 km s . 
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However, more recent studies (Kitayama & Ikeuchi 2000; Dijk- 
stra et al. 2004) find that at high-redshifts (z ^ 3), self-shielding 
and increased cooling efficiency could be strong countering ef- 
fects. These calculations, however, assume spherically symme- 
try, leaving open the possibility of strong feedback for a halo 
with non-isotropic gas profile, illuminated along a low-column 
density line of sight (see Shapiro et al. 2004 for a detailed treat- 
ment of three-dimensional gas dynamics in photo-heated low- 
mass halos). 

Such a drop in the SFR could be detected in the high-redshift 
extension of the 'Lilly - Madau' diagram (Lilly et al. 1996; 
Madau et al. 1996), by directly counting faint galaxies (Barkana 
& Loeb 2000). In practice, the low-mass galaxies susceptible 
to the reionization suppression are faint and suffer from sur- 
face brightness dimming. Whether or not these galaxies will 
fall above the detection threshold of a deep future survey, so 
that they can reveal the effect of reionization, depends crucially 
on the redshift of reionization, the size of the affected galax- 
ies and their typical star-formation efficiencies, as well as the 
amount of dust obscuration. Alternatively, by analyzing the Ly- 
man a absorption spectra of SDSS quasars, Cen & McDonald 
(2002) suggest that we may already have detected a drop in 
the SFR at z ~ 6, through the non-monotonic evolution of the 
mean IGM opacity. In order to improve on this current, low 
signal-to-noise result, deep, high-resolution spectra of bright 
quasars would be required from beyond the epoch of reioniza- 
tion at z > 6 (this may be possible with JWST, see Haiman & 
Loeb 1999). 

Other events that may trace out the early SFR, such as GRBs 
and SNe, could be better suited to detect the reionization drop, 
as they are bright, and could be used as tracers of the SFR in 
arbitrarily faint host galaxies, out to redshifts higher than the 
host galaxies themselves. This is especially useful, since in hi- 
erarchical structure formation scenarios, the characteristic col- 
lapse scale decreases with increasing redshift. As a result, the 
individual galaxies that dominate the ionizing background will 
likely be undetectable if reionization occurred at high redshift 
(z J> 10). Unfortunately, while GRBs are bright, they are rare. 
Based on the expected Swift GRB detection rates (e.g. Figure 
6 in Mesinger et al. 2005), a drop in the GRB rate associated 
with reionization could be detected only at very low, <~ lc, sig- 
nificance, and only if reionization occurred at z re <, 7-10. On 
the other hand, SNe have the benefits of being both very bright 
(compared to galaxies at the very faint end of the luminosity 
function), and occurring much more frequently (compared to 
GRBs). 

In summary, the purpose of this paper is twofold: (i) to con- 
struct expected detection rates of high-z SNe in SNe surveys 
which could be carried out in the future (e.g. with JWST.); (ii) 
to quantify the prospects of detecting a reionization feature from 
such high-redshift SN rates. 

The rest of this paper is organized as follows. In § 2, we 
present our method for predicting the global SFR and SNR, and 
estimating the fraction of SNe detectable in a future survey. In 
§ 3.1, we present our estimates for the SNe detection rates. In 
§ 3.2, we estimate the significance at which the reionization- 
drop in the SNR may be detectable. Finally, we summarize the 
implications of this work and offer our conclusions in § 4. 

Throughout this paper, we assume standard concordance cos- 
mological parameters, (fl\, ^m, fib, n, erg, Ho) = (0.73, 0.27, 
0.044, 1, 0.9, 71 km s" 1 Mpc" 1 ), consistent with the recent mea- 
surements of the power spectrum of CMB temperature anisotropics 



by the Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite (Spergel et al. 2003). 

2. MODELING METHODS 

2.1. The Intrinsic Star Formation and Supernova Rates 

Here we briefly outline our method for obtaining the SFR 
and the intrinsic SNR, using extended Press-Schechter formal- 
ism (EPS) (see Lacey & Cole 1993). The procedure we use 
to obtain the SFR is based on a relatively standard approach 
(see, e.g., Mesinger et al. 2005 and references therein). Readers 
interested in more details are encouraged to consult Mesinger 
et al. (2005). 

We estimate the global SFR density at redshift z as 
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where is the efficiency (by mass) for the conversion of gas 
into stars, MdM(dn(> M, z) / dM) is the mass density contributed 
by halos with total (dark matter + baryonic) masses between M 
and M+dM, t(z) is the age of the universe at redshift z, and P(t) 
is the probability per unit time that new stars form in a mass el- 
ement of age t = f(z)-f(z') (normalized to J Q drP(T) = 1). 
We adopt the fiducial value of e* =0.1 (see, e.g., Cen 2003). 
Note that the star formation efficiency in minihalos (i.e. ha- 
los with virial temperatures below 10 4 K) that contain pristine 
metal-free gas could be significantly lower than 0.1, as sug- 
gested by numerical simulations of the first generation of stars 
(Abel et al. 2002; Bromm et al. 2002). The expected pre- 
reionization SNe rates would then lie closer to our T NU ;> 10 4 K 
curves prior to reionization (see discussion below), making the 
detection of the reionization feature significantly more difficult. 
However, the efficiency is likely to be very sensitive to even 
trace amounts of metalicity (Bromm et al. 2001), and conditions 
for star-formation may result in a standard initial mass function 
(IMF) in gas that has been enriched to metaliticites above a frac- 
tion 10" 4 of the solar value. Indeed, it is unlikely that metal- 
free star-formation in minihalos can produce enough ionizing 
photons to dominate the ionizing background at reionization 
(e.g. Haiman et al. 2000; Haiman & Holder 2003; Sokasian 
et al. 2004). We assume further that star-formation occurs on 
an extended time-scale, corresponding to the dynamical time, 
fdyn ~ [Gp(z)]~ 1/2 (Cen & Ostriker 1992; Gnedin 1996): 
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where p(z) ~ A c p CI j t (z) is the mean mass density interior to 
collapsed spherical halos (e.g. Barkana & Loeb 2001), and 
A c is obtained from the fitting formula in Bryan & Norman 
(1998), with A e = 18tt 2 w 178 in the Einstein-de Sitter model. 
The minimum mass, M m j n (z) in eq. (1), depends on the effi- 
ciency with which gas can cool and collapse into a dark mat- 
ter halo. Prior to reionization and without molecular hydrogen, 
M m i n (z) corresponds to a halo with virial temperature, r v ; r ~ 10 4 
K; with a significant H2 abundance, the threshold decreases to 
T v i r <~ 300 K (Haiman et al. 2000; we use the conversion be- 
tween halo mass and virial temperature as given in Barkana 
& Loeb 2001). Post reionization, the Jeans mass is raised, so 
M m i n (z) could increase. The degree of self-shielding, the abil- 
ity of the halo gas to cool, as well as the amount of H2 present 
in the high-redshift low-mass halos is uncertain (Dijkstra et al. 
2004), and so below we present results for several values of 
M m i n (z), which we will henceforth express in terms of 7^ r (z). 
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Next, from p*(z) we obtain the intrinsic differential SNR 
(number of core collapse SNe per unit redshift per year) with 



dN_ 

dz 



■ »7sn- 



1 dV(z) 
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l+z dz 

where the factor 1/(1 +z) accounts for time dilation, dV(z)/dz 
is the comoving volume in our past light cone per unit redshift, 
and r/sN is the number of SNe per solar mass in stars. For a 
fiducial Salpeter initial mass function (IMF), we obtain 7/ SN ~ 
1/180 Mq 1 , assuming all stars with masses 9 Mq ^ M <> 40 Mq 
become core collapse SNe (Heger et al. 2003). We neglect the 
lifetime of these high mass stars in determining our SNR; this is 
a reasonable assumption as the lifetimes (as well as the spread 
in the lifetimes) are shorter than a unit redshift interval for red- 
shifts of interest. Note that an alternative extreme shape for the 
IMF, consisting entirely of 100-200 Mq stars (Abel et al. 2002; 
Bromm et al. 2002) would yield a similar value for ?7sn- 

We present our SFR densities (top panel) and SNRs (bot- 
tom panel) in Figure 1. The curves correspond to redshift- 
independent virial temperature cutoffs of T vu = 300, 10 4 , 4.5 x 
10 4 , and 1.1 x 10 5 K (or circular velocities of v c j rc = 3, 17, 35, 
and 55 km s _1 , respectively), top to bottom, spanning the ex- 
pected range (Thoul & Weinberg 1996; Dijkstra et al. 2004). 
Also shown are recent results from GOODS (Giavalisco et al. 
2004): the bottom points assume no dust correction and the 
top points are dust corrected according to (Adelberger & Stei- 
del 2000); the statistical error bars lie within the points (Adel- 
berger & Steidel 2000). As there are large uncertainties associ- 
ated with dust correction, each pair of points (top and bottom) 
serves to encompass the expected SFR densities. 

Our SFRs are consistent with other theoretical predictions 
(e.g. Somerville et al. 2001; Barkana & Loeb 2000; Bromm & 
Loeb 2002), as well as other recent estimates from the Hubble 
Ultra Deep Field (Bunker et al. 2004) and the FORS Deep Field 
on the VLT (Gabasch et al. 2004), after they incorporate a factor 
of 5-10 increase in the SFR (Adelberger & Steidel 2000) due to 
dust obscuration. Furthermore, we note that our SNRs, which 
at z <; 5 yield 0.3 - 2 SNe per square arcminute per year, are 
in good agreement with the ~ 1 SN per square arcminute per 
year estimated by Miralda-Escude & Rees (1997) by requiring 
that high-redshift SNe produce a mean metalicity of <~ 0.01 Z 
by z ~ 5. Note, however, that the rates we obtain are signifi- 
cantly higher (by a factor of ~ 60 - 2000 at z <~ 20) than those 
recently found by Wise & Abel (2005). The reason for this 
large difference is that Wise & Abel consider star-formation 
only in minihalos, and they assume a very low star-formation 
efficiency of a single star per minihalo, as may be appropriate 
for star-formation out of pristine (metal-free) gas in the first 
generation of minihalos (Abel et al. 2002; Bromm et al. 2002). 
In contrast, we assume an efficiency of e* = 0.1, which may 
be more appropriate for star-formation in pre-enriched gas that 
dominates the SFR just prior to reionization (including star- 
formation in minihalos). 

As mentioned above, by increasing the cosmological Jeans 
mass, reionization is expected to cause a drop in the SFR (and 
hence the SNR), with the rates going from the horizontally 
striped region in Figure 1 at z > z re to the vertically striped re- 
gion atz < z re . 

The redshift width of this transition is set by a combination 
of large-scale cosmic variance, radiative transfer, and feedback 
effects. For the majority of the paper, we use Az re ~ 1 as a 
rough indicator of the width of the transition we are analyz- 
ing. We distinguish between "reionization" and a "reionization 



feature", and use Az re as an indicator of the width of the later. 
Even with an extended reionization history (Az ~ 10), fairly 
sharp (Az re <, 3) features are likely. We postpone a more de- 
tailed discussion until § 3.2.4. 

The other important factor determining the usefulness of the 
method proposed here is the factor by which the SFR drops 
during the reionization epoch. The size of this drop is medi- 
ated by the effectiveness of self-shielding and gas cooling dur- 
ing photo-heating feedback: i.e. on whether or not the star- 
formation efficiency is significantly suppressed in those halos 
that dominate the SFR and SNR immediately preceding the 
reionization epoch. Given the uncertainties about this feedback 
discussed above, we will consider a range of possibilities be- 
low, parameterized by the modulation in the virial temperature 
threshold for star-formation during reionization. 




FIG. 1. — Upper Panel: SFR densities obtained from our model. The curves 
(top to bottom) correspond to different lower cutoffs on the virial temperatures 
of star-forming halos, T vlr > 300, 10 4 , 4.5 X 10 4 , and 1.1 X 10 5 K (corre- 
sponding to circular velocity thresholds of v c jr C > 3, 17, 35, and 55 km s~'). 
Dots indicate results from GOODS (Giavalisco et al. 2004): the lower set of 
points assume no dust correction, while the upper set of points are dust cor- 
rected; statistical l-cr error bars lie within the points. The lines connecting 
each pair of points span the expected range of SFR densities. Lower Panel: 
the SNR densities accompanying the SFR densities in the top panel. A reion- 
ization feature at redshift z re is expected to be evidenced by a drop in the above 
rates, with rates changing from the horizontally striped region at z > z re to the 
vertically striped region at z < z re , within a transition period lasting for Az re . 



2.2. High Redshift SNe in Future Surveys 

Given the intrinsic star-formation and SN rates, our next task 
is to estimate the number of SNe that could be revealed by a 
future SNe search. In general, the number of SNe per unit red- 
shift, dN exp /dz, that are bright enough to be detectable in an 
exposure of duration f exp can be expressed as 



flWexp dN [ x 
— ; — = —r- I Jsm> hbs,Z) at, 
dz dz Jo 



obs 



(4) 



where (dN / dz)dt b s is the number of SNe which occurred be- 
tween f bs and f bs + c/fobs ago (per unit redshift; note that the 



4 



global mean SNR will evolve only on the Hubble expansion 
time-scale, and can be considered constant over several years), 
and /sn(> f bs,z) is the fraction of SNe which remain visible 
for at least f obs in the observed frame. Then the total number of 
SNe detected in a survey of duration f surv is 

_ f surv Ar; F ov 

surv "2f exp 4ir ^ exp ' Pj 

where Afi FO v is the instrument's field of view, and f SU rv/(2fexp) 
is the number of fields which can be tiled in the survey time, 
fsurv (we add a factor of 1/2 to allow for a second pair of filters 
to aid in the photometric redshift determination; note that this 
provides for imaging in 4 different JWST bands; see discussion 
below). Note also that equation (5) is somewhat idealized, in 
that it assumes continuous integration for the duration of the 
survey, and e.g., does not account for time required to slew the 
instrument to observe different fields. In principle, each field 
has to have repeated observations (to detect SNe by their vari- 
ability), and therefore any dedicated survey should target fields 
that have already been observed. Furthermore, a long, dedi- 
cated program may not be necessary, because the effect could 
be detected with relatively few fields (at least under optimistic 
assumptions; see discussion below), and several fields with re- 
peated imaging (separated by > 1-2 years) may already be 
available from other projects; these fields can then be used for 
the SN search. 

Additionally, it might prove useful to consider the luminosity 
of a SN host galaxy as a means of biasing the sample towards 
the small mass halos which are most susceptible to the reion- 
ization suppression. Ignoring SNe which originate in galaxies 
above a certain threshold luminosity at a given redshift could 
limit the sample to those small halos which are most affected 
by reionization, and hence decrease the "noise" component in 
the reionization signal (e.g. eq. (6)). If there turns out to be 
strong correlation between SN type and host halo properties 
(e.g. due to metalicity), then, in principle, such correlations 
can be empirically discovered and used to further narrow the 
SN sample to those most susceptible to feedback - even if the 
expected correlation cannot be predicted ab-initio. 

In general, /sn(> ?obs,z) in equation (4), i.e. the fraction of 
SNe which remain visible for at least f Q b s , depends on (i) the 
properties of the SN, in particular their peak magnitude and 
lightcurve, and the distribution of these properties among SNe, 
and (ii) on the properties of the telescope, such as sensitivity, 
spectral coverage, and field of view. In the next two subsec- 
tions, we discuss our assumptions and modeling of both of these 
in turn. 

2.2.1. Empirical Calibration of SN Properties 

At each redshift, we run Monte-Carlo simulations to deter- 
mine /sn(> fobs,z) in equation (4). We use the observed prop- 
erties of local core-collapse SNe (CCSNe) in estimating /sn(> 
tobs,z). For the high redshifts of interest here, we only consider 
core collapse SNe of Type II. SNe resulting from the collapse 
of Chandreshekar-mass white dwarfs (Type la) are expected 
to be extremely rare at high redshifts (z <; 6), as the delay be- 
tween the formation of the progenitor and the SN event ( ^ 1 
Gyr; Strolger et al. 2004) is longer than the age of the universe 
at these redshifts. Local core CCSNe come in two important 
varieties, types IIP and IIL, differentiated by their lightcurve 
shapes. We ignore the extremely rare additional CCSN types, 
e.g Type Iln and lib, which appear to have significant interac- 
tion with circumstellar material and constitute less than 10% of 



all CCSNe. Type Ib/c, which may or may not also result from 
core collapse, have luminosities and light-curves that are sim- 
ilar to Type IIL and occur less frequently. While the relative 
numbers of Type IIP and Type IIL SNe are not known even for 
nearby SNe, recent estimates imply that they are approximately 
equal in frequency (Cappellaro et al. 1997). We therefore as- 
sume that 50% of the high-redshift SNe are Type IIP and 50% 
are Type IIL. 

CCSNe result from the collapse of the degenerate cores of 
high-mass stars. The luminosity of CCSNe is derived from the 
initial shock caused by the core-collapse which ionizes material 
and fuses unstable metal isotopes (see Leibundgut & Suntzeff 
(2003) and references therein for a more detailed description 
of SN lightcurves). In the early stages of the SN, the shock 
caused by the core collapse breaks out from the surface of the 
progenitor (typically high mass red giants), resulting in a bright 
initial peak in the light curve that lasts less than a few days 
in the rest-frame of the SN. As the shock front cools, the SN 
dims. However, the SN may then reach a plateau of constant 
luminosity in the light curve, believed to be caused by a wave of 
recombining material (ionized in the shock) receding through 
the envelope. The duration and strength of this plateau depends 
on the depth and mass of the progenitor envelope, as well as the 
explosion energy, with those SNe exhibiting a strong plateau 
classified as Type IIP. A typical plateau duration is <; 100 rest- 
frame days (Patat et al. 1994). In Type IIL SNe, this plateau 
is nearly non-existent, and the lightcurve smoothly transitions 
from the rapid decline of the cooling shock to a slower decline 
where the luminosity is powered by the radioactive decay of 
metals in the SN nebula. After the plateau, Type IIP SNe also 
enter this slowly declining 'nebular' phase. 

These observationally determined behaviors have been sum- 
marized in a useful form as lightcurve templates in Doggett & 
Branch (1985). We use these template lightcurves in determin- 
ing /sn(> t b>s,z), and normalize the lightcurves using gaussian- 
distributed peak magnitudes (i.e. log-normally distributed in 
peak flux) determined by Richardson et al. (2002) from a large 
sample of local Type IIP and Type IIL SNe (see Table 1). We 
perform the Monte-Carlo simulations with both the dust cor- 
rected, and dust uncorrected values in Richardson et al. (2002), 
since the dust production history of the early universe is poorly 
understood and is essentially unconstrained empirically. 

We use a combined high-resolution HST STIS + ground- 
based spectrum of the Type IIP supernova SN1999em (the 'Novem- 
ber 5th' spectrum of Baron et al. 2000) as the template SN 
spectrum in order to obtain TT-corrections (with the Doggett & 
Branch 1985 lightcurves given in the restframe B filter). This 
spectrum was obtained within 10 days of maximum light, i.e. 
during the initial decline of the SN brightness, and has been 
dereddened by Ay = 0.3 mag (Baron et al. 2000; Hamuy et al. 
2001). While the spectrum, and hence the Zf-corrections, of 
SNe evolve during the lightcurve, this effect is not strong for the 
wavelengths of interest (Patat et al. 1994), especially since the 
lightcurve template we use is well matched to the wavelengths 
being probed by the observations we consider below (leading 
to small -corrections). We have also used this Type IIP SN 
template spectrum to calculate TT-corrections for Type IIL SNe. 
This is necessary due to the lack of restframe UV spectra of 
Type IIL SNe that can be combined with optical spectra, and 
justifiable because the /f-corrections are relatively small and the 
broadband colors of both Type IIP and IIL SNe (a measure of 
the spectral shape that determines the Zf-corrections) are sim- 
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ilar, at least in the optical (Patat et al. 1994). Note that we 
assume that very-high redshift core-collapse SNe are similar to 
local SNe in their spectra, peak luminosities, and temporal evo- 
lution. However, these assumption do not significantly impact 
our conclusions below, as long as (i) the average properties of 
the SNe do not change rapidly at high redshift (which could 
mimic the reionization drop), (ii) they do not become preferen- 
tially underluminous (which would make high-z SNe less de- 
tectable, lowering the statistical confidence at which the reion- 
ization drop is measured), and (iii) the detection efficiency does 
not change rapidly with redshift (e.g. due to instrument param- 
eters or spectral lines). 

Table 1 

Means and standard deviations of the adopted peak 
absolute magnitudes. values are taken from 
Richardson et al. (2002). Note: an M b = -17 SN would be 

DETECTABLE OUT TO z ~ 8.2 AT THE FLUX THRESHOLD OF 3 nJy 
AT 4.5 /an. 

Corrected for Dust Not Corrected for Dust 
SNType (M B ) a (M B ) a 

Hp a7W TTT2 AHTl L23 
IIL -18.03 0.9 -17.80 0.88 



o 




-0.5 L_. , , , i , , , , i , , , , i , , , , i , . , , 

-0.5 0.0 0.5 1.0 1.5 2.0 

3.5 micron - 4.5 micron 



FIG. 2. — Infrared colors of a Type II SN as a function of its redshift. Solid 
curve: The template spectrum of the text, from the Type IIP SN 1999em. Di- 
amonds are placed at intervals of Az = 1, labeled with the redshift. Error bars 
of 0.15 mag are shown to give a sense of the photometric errors that may be 
expected, excluding the many possible systematic effects. Dashed curve: The 
colors of a blackbody at different redshifts (temperatures). The asterisk marks 
the color of a 9000K blackbody at z = 10, triangles are at intervals of Az = 1 
(or at a constant redshift but with a corresponding different temperature) 



In order to use the method outlined above to probe reion- 
ization, the SN redshifts must also be known to an accuracy 
of Az <^ 1. SN redshifts can be determined via spectroscopy 
of either the SN itself, or of the host galaxy. However, as we 
have already noted, the host galaxies may only be marginally 
detectable even in imaging, and the SNe may be too faint for 
anything other than extremely low resolution spectroscopy. We 
present here only a very brief example of the possibility of ob- 
taining redshifts from the extremely low resolution (A/AA ~ 5) 
spectra provided by multiband imaging: a complete investiga- 
tion of this possibility is warranted, but is beyond the scope of 
this paper. 

To the extent that Type II SNe spectra can be represented as a 
sequence of blackbodies of different temperatures (e.g. Dahlen 
& Fransson 1999) photometric redshifts will be impossible to 
obtain without information about the SN epoch, since temper- 
ature and redshift would be degenerate. However, local Type 
IIP SNe show significant deviations from a blackbody in the 
UV (A < 3500A) due to metal-line blanketing in the SN pho- 
tosphere, providing spectral signatures that could be used as 
redshift indicators, depending on their strength. In Figure 2, we 
show the evolution with redshift of the infrared colors (in bands 
accessible with JWST; see below) of our template SN spectrum, 
compared with the color evolution of a blackbody. The figure 
shows that the template spectrum deviates significantly from a 
blackbody. If the spectrum of the SN is always the same as the 
template spectrum, then there are good prospects for obtain- 
ing photometric redshifts for these SNe, at least in the redshift 
range z = 7 — 13. While the figure shows that multiple redshifts 
may be possible at fixed observed colors, the degenerate so- 
lutions would correspond to z > 16 SNe; contamination from 
such high redshift will be mitigated by the fact that these SNe 
are likely to be too faint to be detected. Of course, more de- 
tailed studies of the UV behavior of local SNe, especially their 
variety and spectral evolution, will be necessary to confirm the 
possible use of photometric redshifts. 



2.2.2. SNe Detectability and Survey Parameters 

As a specific example for the number of detectable SNe in a 
future sample, we consider observations by JWST, a 6m diame- 
ter space telescope, scheduled for launch in 2013. Similar con- 
straints could be delivered by the Joint, Efficient Dark-energy 
Investigation (JEDI) 1 (Wang et al. 2004), a proposed instru- 
ment which will require a <~ 100 times longer integration time 
to reach the same detection threshold, but has a ~ 100 times 
larger FOV. The relevant instrument on JWST is NIRcam, 2 a 
near-infrared imaging detector with a FOV of 2.3' x 4.6'. A 
field can be observed in two filters simultaneously. NIRcam 
will have five broadband filters (with resolution A/AA ~ 5). 
We model the filter response as tophat functions with central 
wavelengths of 1.5, 2.0, 2.7, 3.5, and 4.5 /mi. For concreteness, 
below we will present results only for the 4.5 and 3.5 /im fil- 
ters, since they are the longest-wavelength JWST bands; how- 
ever, we allow time for imaging in two other bands, if needed 
for photometric redshift determinations. The current estimate 
of the JWST detection threshold at 4.5 /xm is Z 3 nJy for a 10 
a detection and an exposure time of 10 5 s. The 3.5 /mi band is 
more sensitive, with a detection threshold of > 1 nJy for a 10 
a detection and an exposure time of 10 5 s. 

We show our results for /sn(> ?obs,z) in Figure 3. The solid 
curves correspond to z = 7, the dashed curves correspond to 
z = 10, and the dotted curves correspond to z = 13. In each 
panel, the top set of curves assumes a flux threshold of 3 nJy 
(or an exposure time of f exp = 10 5 seconds in the 4.5 /mi band), 
and the bottom set assumes 9.5 nJy (r exp = 10 4 seconds in the 
4.5 /im band, background dominated). The top panel further 
assumes no dust extinction, while in the bottom panel, we adopt 
the same dust extinction as in the low redshift sample (Richard- 
son et al. 2002). Understandably, the distributions get wider 
as redshift increases (due to time dilation), but the total visi- 
ble fraction, /sn(> 0,z), gets smaller (due to the increase in 

'http://jedi.nhn.ou.edu/ 

2 See http://ircamera.as.arizona.edu/nircam for further details. 
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luminosity distance). The double bump feature in some of the 
curves corresponds to the plateau of Type IIP SNe lightcurves 
discussed above. 




5xl0 3 



1.5xl0 3 



2xl0 3 
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(days) 



FIG. 3. — Fraction of SNe which remains visible for an observed duration 
of ; bs or longer. The curves correspond to SN redshifts of z = 7 (solid curve), 
z = 10 (dashed curve), and z = 13 (dotted curve). In each panel, the top set 
of curves assumes a flux threshold of 3 nJy (exposure time of t ex p = 10 5 s with 
the 4.5 fim JWST), and the bottom set assumes 9.5 nJy (? e xp = 10 4 s). The top 
panel assumes no dust extinction, while the bottom panel assumes the same 
dust extinction as observed within the low redshift sample (Richardson et al. 
2002). 

As seen in Figure 3, most of the supernovae that are bright 
enough to be visible at all, will remain visible for up to <~ 1 
- 2 years. Hence, in order to catch the most SNe, it will be 
necessary to have repeat observations of the SN survey fields 
a few years apart, to insure that most of the observed SNe will 
be identified as new sources or sources that have disappeared. 
A time between observations that is comparable to or larger 
than the SN duration will be the optimal strategy for detecting 
the most SNe (Nemiroff 2003). Note that we ignore the time 
required to collect reference images, since observations con- 
ducted for other programs will likely provide a sufficient set of 
such reference images. However, as shown in equation (5), we 
allow time for the field to be imaged in four different bands, to 
aid in photometric redshift determination. 

To be more explicit, we find that in order to obtain the largest 
number of high-redshift SNe, in general it is a more efficient 
use of JWSTs time to 'tile' multiple fields rather than 'stare' 
for extended periods of time ( ^> 10 4 s) at the same field (Ne- 
miroff 2003). This is because of strong time dilation at these 
redshifts. As can be seen from Figure 3, most detectable SNe 
will remain above the detection threshold for several months, 
even assuming 10 4 s integration times. Also evident from Fig- 
ure 3 is that the increase in the total visible fraction of SNe 
going from f exp = 10 4 s to f exp = 10 5 s, is less than the factor of 
10 increase in exposure time. As a result, a fiducial 1-yr JWST 
survey would therefore detect more SNe using f exp = 10 4 s than 
using f exp = 10 5 s (see Figures 4a and 5a). Understandably, this 
conclusion does not hold for very high-redshifts, z ^ 14, where 



SNe are extremely faint, and require very long exposure times 
to be detectable. However, even with such long exposure times, 
very few SNe will be detectable at these large redshifts, render- 
ing the use of longer exposure times unnecessary. 

3. RESULTS AND DISCUSSION 

3.1. SNe Detection Rates 

The number of SNe that could be detectable in putative fu- 
ture surveys are shown in Figures 4a and 5a. The curves cor- 
respond to the same virial temperature cutoffs for star-forming 
halos as in Figure 1. Solid lines assume no dust obscuration; 
dashed lines include a correction for dust obscuration as dis- 
cussed above. Figure 4a shows results assuming flux density 
thresholds of 9.5 nJy (or f exp = 10 4 s with the 4.5 [im JWST fil- 
ter) (top panel) and 3 nJy (/ exp = 10 5 s) (bottom panel). Figure 
5a shows results with the 3.5 /im filter assuming equivalent ex- 
posure times: flux density thresholds of 3.2 nJy (f exp = 10 4 s 
with the 3.5 JWST filter) (top panel) and 1 nJy (f exp = 10 5 s) 
(bottom panel). The right vertical axis displays the number of 
SNe per unit redshift per FOV (2.3' x 4.6'); the left vertical axis 
shows the number of SNe per unit redshift in a fiducial 1-year 
survey. As mentioned above, reionization should be marked by 
a transition from the region bounded by the top two solid curves 
to the region enclosed by the bottom three solid curves (or the 
analog with the dashed curves if dust is present at the time of 
reionization). 

We note that our expected rates are somewhat higher than 
those in Dahlen & Fransson (1999), a previous study which 
included SNe lightcurves and spectra in the analysis. For ex- 
ample, we find 4-24 SNe per field at z ^ 5 in the 4.5 /im filter 
with ? exp = 10 5 s, compared to ~ 0.7 SNe per field at z ^ 5 ob- 
tained by Dahlen & Fransson (1999) (after updating their JWST 
specifications to the current version). However, they use SFRs 
extrapolated from the low-redshift data available at the time, 
which are not a good fit to recent high-z SFR estimates (Gi- 
avalisco et al. 2004; Gabasch et al. 2004; Bunker et al. 2004), 
and are lower than our z Z 5 SFRs by a factor of 6 - 40. 3 Taking 
this factor into account, their procedure yields 4-27 SNe per 
field at z ^ 5, which is in excellent agreement with our estimate 
of 4 - 24 SNe per field at z ^ 5. 

3.2. Detecting Reionization Features 

To quantify the feasibility of detecting the reionization fea- 
ture, we assume a Poisson S/N, 

(A^pre) " (A/post) 



S/N= V«field- 



(6) 



V / (A'pre) + (A r post) ' 

where rcgeid is the number of JWST fields, and (Ap re ) and (Af post ) 
are the number of observable SNe per field per unit redshift 
pre and post reionization, averaged over z re < z < z re + Az re and 
z re - Az re < z < z re , respectively. In our fiducial model we take 
Az r e = 1, but below we also investigate the impacts of a larger 
value of Az re . Equation (6) yields a simple figure of merit that 
does not distinguish between the theoretical means of the distri- 
butions we generate and the future observational samples which 
will have to be interpreted as numbers drawn from distributions 
with unknown means. Likewise, simple ^/N errors are under- 
estimates for small A'. However, we have verified that at our 

3 The possibility that the SFR in a "Lilly-Madau diagram" remains flat, or 
even increases, at redshifts z > 5, owing to star-formation in early, low-mass 
halos, is also expected theoretically (see, e.g., Fig.l in Bromm & Loeb 2002 
and associated discussion). 
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FIG. 4. — (a) Number of high-redshift SNe detectable with the 4.5 /im JWST filter. The curves correspond to the same virial temperature cutoffs as in Figure 1. 
Solid curves assume no dust obscuration; dashed curves adopt dust obscuration in the same amount as observed in the low redshift SNe sample. The figure shows 
results assuming flux density thresholds of 9.5 nJy (or f e x P = 10 4 s with JWST) (top panel) and 3 nJy (or fe xp = 10 5 s with JWST) (bottom panel). The right vertical 
axis displays the number of SNe per unit redshift per field; the left vertical axis shows the number of SNe per unit redshift found in r sll rv /(2fcx P ) such fields (i.e. the 
differential version of eq. (5) with f surv = 1 yr). Reionization should be marked by a transition from the region bounded by the top two solid curves to the region 
enclosed by the bottom three solid curves (or the analog with the dashed curves if dust is present at the time of reionization). (b) Number of 2.3' X 4.6' JWST fields 
required to detect a reionization feature occurring at z re , with a S/N > 3. Solid curves correspond to no dust obscuration; dashed curves include dust. The top 



panels assume ?exp = 10 4 s; bottom panels assume ; e xp = 10 5 s. Left: Reionization drops in the SNR analogous to T vu > 300 K — * 10 4 K, r v ; 
K, and T vk > 300 K -* 1.1 X 10 5 K, (top to bottom). Right: Reionization drops in the SNR analogous to 7 vlr > 10 4 K -> 4.5 X 10 4 K, and T v 
K, (top to bottom). 



> 300 K - 
• > 10 4 K - 



• 4.5 X 10 4 
» 1.1 x 10 5 




FIG. 5. — Same as Figure 4, but with the 3.5 filter, instead of the 4.5 (im one. Note that we present results for comparable exposure times, hence the sensitivity 
thresholds are three times lower than in Figure 4, due to the disparate sensitivities of the 3.5 (im and 4.5 (im filters. 
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chosen limiting value of S/N <; 3, these issues have little im- 
pact on the results. A full statistical treatment (e.g. comparing 
redshift distribution of counts with a Fisher matrix technique; 
Holder et al. 2001) is beyond the scope of this paper, and would 
be of little value given the uncertainties in our estimates. How- 
ever, we note that the statistics could be improved by combin- 
ing several (e.g. Az = 0.5 wide) bins prior to and following 
reionization, or by applying the null-hypothesis of a smoothly 
evolving SNR directly to the unbinned table of the observed 
SNe (Holder et al. 2001). 

In Figures 4b and 5b, we plot the number of 4.5 and 3.5 
/im JWST fields, respectively, required to detect the reionization 
feature at z re with a S/N <; 3. Solid curves correspond to no 
dust obscuration; dashed curves include dust corrections. The 
top panels assume f exp = 10 4 s; bottom panels assume f exp = 10 5 
s. The left panels show the results assuming that reionization 
modifies the virial temperature cutoff as r v ; r ^ 300 K — > 10 4 
K, T vir ^ 300 K -> 4.5 x 10 4 K, and T vir ^ 300 K -> 1 . 1 x 10 5 
K, from top to bottom. The right panels show results under the 
less pronounced transitions of T v \ r <; 10 4 K — ■> 4.5 x 10 4 K, and 
r vir £ 10 4 K -> 1 . 1 x 10 5 K, from top to bottom. To obtain the 
corresponding number of necessary SNe, one can multiply the 
number of required fields in Figures 4b and 5b with the corre- 
sponding number of SNe per unit redshift per field in Figures 4a 
and 5a, respectively. Typically, the detection of the reionization 
drop with S/N ^ 3 requires at least (N pK ) ^ 10 SNe. 

The relative flatness of most of the curves in Figures 4b and 
5b over the ranges 7 <, z re 13 and 6 <, z re ^ 10, respectively, 
is due to the fact that even though the number of detections de- 
creases with increasing redshift (making the detection of the 
reionization drop more difficult), the separation between the 
various suppression curves increases (making the detection of 
the reionization drop easier). We note from the figures that 
detecting the reionization drop in the SNe distributions con- 
structed from the 3.5 /zm filter is more efficient (requires fewer 
comparable fields) than using the 4.5 pm filter at redshifts 6 ^ 
z re ^ 10. At redshifts larger than z re ~ 10, the 4.5 /im filter be- 
comes more efficient, and can be used to detect the reionization 
drop out to z re <S 15 in a 1 year survey, in the most optimistic 
scenario. However, such conclusions about SNe at redshifts as 
high as z re <; 13-15 are highly uncertain: (i) there could be 
an intrinsic physical cut-off to the bright-end tail of the SNe 
luminosity distributions, drastically reducing the number of de- 
tectable SNe at very high redshifts, and/or (ii) detections in the 
3.5 /im band, in addition to the the 4.5 /im band, could be neces- 
sary for photometric redshift determination (see § 2.2. 1), which 
could further reduce the number of very high redshift SNe with 
accurate redshift determinations. 

In the most optimistic scenario we considered (with the thresh- 
old for star formation raised from r v ; r (z > z re ) ^ 300 K to r vrr (z < 
Zre) <; 1.1 x 10 5 K), corresponding to maximum photo-heating 
feedback (abundant H2 cooling and no self-shielding), the drop 
can be detected with just 8 (3.5 /im) fields over the range 6 ^ 
Zre & 10 and 20 (4.5 /im) fields over the range 6 ^ z re ^ 13 with 
fexp = 10 5 s. Using shorter exposures, it can be detected with 20 
(3.5 /im) fields over the range 6 <> z re < 10 and 100-300 (4.5 
/im) fields over the range 6 <J z re ^ 13 with f exp = 10 4 s. Less 
optimistic scenarios can still allow for a detection of the reion- 
ization drop with only tens of fields. However, a reionization 
feature at z re 15, would require ^ 1000 fields in order for the 
associated drop to be detected. 

In the most pessimistic scenario (with the threshold changing 



from r v ; r <; 10 4 K — > 10 4 K), corresponding to very strong self- 
shielding and no H2 at reionization (i.e. the complete absence of 
any photo-heating feedback), reionization would not affect the 
SNR and could not be detected. Nevertheless, if a sharp reion- 
ization feature is detected independently (through Lya absorp- 
tion spectra of galaxies or quasars, 21 -cm signatures, or CMB 
anisotropies), the lack of a detected drop in the SNR would pro- 
vide valuable evidence that low-mass halos can withstand the 
effects of photo-ionization heating. Note that without a compli- 
menting detection of a sharp reionization feature, one can not 
distinguish between a very smooth, gradual reionization and an 
absence of photo-heating feedback. 

Finally, we comment on the origin of the spike at z re ~ 6 
in the curves in Figure 4b. Apparently, the detection of the 
reionization feature in the 4.5 /mi band at z re ~ 6 - 7 has an 
additional difficulty. This is caused by a combination of two 
effects: (i) the decreasing separation between the observable 
SNe distributions (as seen in Figure 4a) at lower redshifts, and 
(ii) a large-equivalent-width Ha emission line in the spectrum 
of our template SNe, SN1999em, at the rest frame wavelength 
X Ha rs 6563 A. (The latter effect also produces the slight bump 
in the curves in Figure 4a at the redshift in which the 4.5 /im fil- 
ter probes the line, zhq 4.5^m ~ 4.5/im/AH Q - 1 ~ 6.) The 
strong Ha emission line influences the TT-correction so as to in- 
crease the number of detectable SNe at z Hq 4.5 Mm - This reduces 
the apparent drop in the SFR, if the reionization feature oc- 
curs at redshifts just above zna 4.5/^m (and, likewise, increases 
the size of the drop if the reionization feature occurs just be- 
low zhc 4.5 J um)- The detectability of the reionization feature is 
only affected at z Hq 4.5 Mm ~ 6 due to the combination of (i) and 
(ii) above. The presence of the Ha line is nevertheless po- 
tentially important, because reionization could have occurred 
around this redshift (Mesinger & Haiman 2004). This high- 
lights the need for an improved statistical description of SNe 
template spectra from low-redshift data, allowing one to "model 
out" an offending emission line, or for observations of the high- 
z SNe in a different filter (as can be seen from Figure 5, using 
the 3.5 /im filter bypasses this issue). 

In order to obtain high-redshift SNRs, several simplifica- 
tions have been made in our analysis, involving the poorly con- 
strained redshift evolution (more precisely, the lack of evolu- 
tion) of the stellar initial mass function (IMF), SN progenitors, 
star formation efficiency, SN lightcurves and spectra, etc. Be- 
cause of this, the numbers presented in Figures 4 and 5 should 
be regarded only as rough estimates of what a future survey 
might achieve. However, our main conclusions rely on a sud- 
den feature in the redshift evolution of the SNR, and should 
remain valid, as long as evolution in other parameters occurs 
over a much wider redshift range than the reionization feature 
(Wyithe & Loeb 2004a; Barkana & Loeb 2004b), and as long as 
the pre-reionization SNR is not significantly reduced (since the 
latter would diminish the confidence at which any subsequent 
drop in the SNR, caused by reionization, could be detected). 

3.2.1. Pop-Ill SNe 

As discussed above, high-redshift SNe whose progenitor stars 
are formed from metal-free gas within minihalos could be in- 
trinsically very different from the low-redshift SNe, due to dif- 
ferences in the progenitor environments (e.g. very low met- 
alicities; Abel et al. 2002; Bromm et al. 2002). If such dif- 
ferences could be identified and detected, then these "pop-Ill" 
SNe could provide valuable information about primordial stars 
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and their environments. Indeed, Wise & Abel (2005) recently 
studied the redshift-distribution of such primordial SNe, but 
only briefly addressed the issue of their detectability. In order 
to directly assess the number of such SNe among the hypo- 
thetical SNe samples we obtained here, in Figure 6, we plot the 
fraction of SNe whose progenitor stars are located in minihalos. 
The solid curve assumes our fiducial model with a Salpeter IMF 
and e* m i n ihaio = 0. 1 ; the dotted curve assumes that each minihalo 
produces only a single star, and hence a single SN, over a dy- 
namical time (assuming that strong feedback from this star dis- 
rupts any future star formation; as in Wise & Abel 2005). The 
figure shows that with an unevolving star formation efficiency, 
progenitor stars in minihalos would account for over half of the 
SNe at z J; 9. On the other hand, in the extreme case of a single 
SNe per minihalo, progenitor stars in minihalos would account 
for less than ~ 1% of the z ~ 10 SNe (although at the earliest 
epochs, z ^ 22, they would still constitute over half of all SNe). 
Given the overall detection rate of several hundred z k5 10 SNe 
with a 1 year JWST survey we have found above, even in this 
extreme case, Figure 6 implies that several of these SNe could 
be caused by pop-Ill stars. 
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FIG. 6. — Fraction of SNe whose progenitor stars are located in miniha- 
los. The solid curve assumes our fiducial model with a Salpeter IMF and 
e*minihalo = 0.1; the dotted curve assumes that each minihalo produces only a 
single star (and therefore a single SN) over a dynamical time. 



3.2.2. Varying Minihalo Star-Formation Efficiency 

The uncertainties in the typical star-formation efficiency in 
minihalos can also directly influence our estimates for the de- 
tectability of the reionization feature. We have already shown 
that the feature is detectable from SNe distributions even in the 
absence of star formation in minihalos, but only provided that 
halos do not self-shield effectively and the virial temperature 
for star-formation is raised well above 10 4 K (e.g. see the right 
panels of figures 4b and 5b). If self-shielding is very effective, 
the reionization drop is determined by the evolution of the SFR 
in minihalos. In Figure 7, we plot the S/N with which the reion- 
ization feature is detectable from a fiducial one-year survey, as 
a function of the assumed (constant) star-formation efficiency 
in minihalos, e* m i n jh a io- The three curves correspond to z re = 7 
(solid curve), 10 (dotted curve), 13 (dashed curve). All curves 
assume the transition r vir (z > z re ) ^ 300 K — > r vir (z < z re ) ^ 10 4 
K, as well as a flux density threshold of 3 nJy (f exp = 10 5 s with 
the 4.5 /im JWST filter), and no dust extinction. The figure 
shows that in this case, the reionization signal is detectable at 
6 ^ z re ^ 13 with S/N > 2 as long as the star formation effi- 
ciency in minihalos exceeds e* m i n jh a io ~ 0.05. 
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FIG. 7. — S/N with which reionization can be detected as a function of 
the assumed (constant) star— formation efficiency in minihalos, e* m inihalo- 

The 

three curves correspond to three different reionization redshifts z re = 7 (solid 
curve), 10 (dotted curve), and 13 (dashed curve). All curves assume the tran- 
sition r vir (z > z rc ) > 300 K — > r vir (z < z re ) Z !0 4 K, as well as a flux density 
threshold of 3 nJy (< exp = 10 5 s with the 4.5 fim JWST filter), no dust extinction, 
and a « surv = lyr survey (or, correspondingly, using lyr/(2« eX p) images with a 
2.3' X 4.6' FOV). 



3.2.3. Varying Az re 

As mentioned previously, the width of a reionization feature, 
which we have parameterized by Az re (see the clarification be- 
low eq.(6)), is uncertain. In our fiducial model we have as- 
sumed to study fairly sharp features with Az re = 1- Here we 
investigate the detectability of more extended reionization fea- 
tures. In Figure 8, we plot the maximum reionization redshift, 
Zre max, up to which the reionization signal would be detectable 
at S/N <; 3 in a 1 yr survey, as a function of Az re • To approx- 
imate the redshift-smearing due to any self-regulation or scat- 
ter, we simply follow the procedure discussed above (follow- 
ing eq.6), but take averages of the SNR in increasingly wider 
redshift bins prior to and following z re . This effectively aver- 
ages the SNR in two bins whose mean redshifts are separated 
by Az re and mimics a suppression that is spread-out over this 
redshift, rather than instantaneous. A reionization feature oc- 
curring in the range 6 z re Ss z re max, smeared over a redshift 
range of Az re , would be detectable out to the redshifts shown in 




FIG. 8. — Maximum reionization redshift, z re raa x, up to which the reioniza- 
tion signal would be detectable at S/N > 3 in a 1 yr survey, as a function of 
Az re . The curves correspond to the reionization transitions of T v ; r > 300 K — > 
10 4 K, r vir > 300 K -> 4.5 X 10 4 K, and T vlr > 300 K -> 1.1 X 10 5 K, (left to 
right). Solid curves correspond to no dust obscuration; dashed curves include 
dust. All curves assume a detection threshold of 1 nJy, analogous to fc x p = 10 5 
s exposures in the 3.5 ^m /WSJ band. 
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Figure 8. The curves correspond to the reionization transitions 
of T vir £ 300 K -» 10 4 K, T vir £ 300 K -> 4.5 x 10 4 K, and 
T vir £ 300 K -» 1.1 x 10 5 K, (left to right). Solid curves cor- 
respond to no dust obscuration; dashed curves include dust. 
All curves assume a detection threshold of 1 nJy, analogous to 
f eX p = 10 5 s exposures in the 3.5 JWST band. From the fig- 
ure, one can note that the most optimistic scenario (T vlv <; 300 
K — > 1.1 x 10 5 K), a reionization feature is detectable even if 
it occurs over a redshift range Az re ~ 4. In other cases, the 
detection of the reionization feature requires Az re & 1-3. 

3.2.4. Feasibility of a Fairly Sudden Drop in SNe Rates 

The details and duration of the reionization epoch, and hence 
the shape and width of the drop in SNe rates, are unknown at 
this time. Results from WMAP and the SDSS QSOs offer sug- 
gestive, albeit not conclusive, evidence that reionization was 
extended in redshift. 4 If reionization is indeed quite extended 
in redshift (Az <~ 10), the associated drop in the SFR and SNR 
could be smeared out too much to be detectable. However, even 
in the pessimistic (as it pertains to our analysis below) scenario 
where some effective width of the reionization feature is as 
large as Az <~ 10, the shape of the feature need not be smooth, 
and could contain "sharp" (Az re <~ 1) drops. Such transitions 
are probable in reionization histories in which different sources 
dominate different epochs (e.g. Haiman & Holder 2003). 

The relevant process in determining the width and shape of 
the drop in SNRs is the evolution of the volume filling factor 
of ionized regions and their correlation with the small, vulner- 
able halos 5 whose SFRs are sensitive to the thermal state of the 
IGM. In other words, the sharpness of the drop in SNRs, in each 
epoch during reionization, will depend on: 

1 . the nature of the dominant ionizing sources 

2. the ionizing efficiency of the dominant ionizing sources 

3 . the level of synchronization of small halo formation with 
the formation of dominant ionizing sources 

4. the level of synchronization of the dominant ionizing 
sources 

If the small, vulnerable halos themselves are the dominant 
ionizing sources at that particular epoch, then a sharp reioniza- 
tion feature could result if: (2) is high, and (3) is moderate (in 
this case (3) and (4) are the same). If (3) is too high (i.e. the 
formation of small, vulnerable halos is very clustered in time 
and space), then negative feedback from the ionizing radiation 
could delay substantial growth of HII bubbles until small halos 
are no longer forming prodigiously enough to serve as signposts 
for reionization. If (3) is too low, there might not be enough af- 
fected halos to notice the suppression; or large, isolated patches 
might have to wait a long time for their own ionizing sources 
to form (if (2) is not very high), thus smearing out the signal 
through pure cosmic variance. 

4 Note that there are two, often confused processes and time-scales associated 
with reionization: (i) the increase in intensity of the ionizing background and 
the mean free path of ionizing photons and (ii) the increase in the filling factor 
of ionized regions. For the purposes of this paper, we concern ourselves with 
(ii). 

5 Note the distinction here between "small halos" and "minihalos". In the dis- 
cussion below, we use the term small halos to denote all halos whose SFRs will 
be suppressed by reionization. Thus small halos include minihalos as well as 
larger halos, depending on their susceptibility to negative photo-heating feed- 
back. 



If the small, vulnerable halos are not the dominant ionizing 
sources (as would be expected for the later periods of an ex- 
tended reionization), then a sharp feature could result if (2) is 
high enough to reasonably counter cosmic variance. However, 
we could relax our fine-tuning on (3) above, since feedback no- 
longer hinders the growth of HII regions. There merely need 
to be enough small halos at that epoch to act as signposts for 
reionization. From Figure 6, we see that this a reasonable as- 
sumption, especially given the fact that most small halos which 
are still forming at such late stages are probably not going to be 
very near the large overdensities which were likely to be ion- 
ized during earlier stages (Furlanetto & Oh 2005; Ricotti et al. 
2002). We also require (4) to be reasonably high (i.e. that the 
dominant ionizing sources appear around the same time, with- 
out too much cosmic scatter). Below, we further quantify such 
a scenario. 

One can get a sense of the possible shapes of the reionization 
feature through an estimate of the evolution of the filling factor 
of ionized regions, Fnn(z), (c.f. Barkana & Loeb 2001; Haiman 
& Holder 2003): 



dF ml (z) 
dt 



■■ e*/es 



^Vph/b^coi(>M min (z),z) 



Vph/b 

0.76 



dt 



-a B C(n Q H )(l+z) 3 F ml . 

(7) 

Here / esc is the escape fraction of ionizing photons, Af ph / b is 
the number of ionizing photons per baryon emitted by a typical 
source, F co i(> M,z) is the fraction of baryons that reside in col- 
lapsed halos with a total mass greater than M at redshift z, as is 
the hydrogen case B recombination coefficient, C = (n^)/(n H ) 2 
is the clumping factor, and (n^) is the current hydrogen number 
density. The first term on the right hand side accounts for "new" 
ionizations contributing to the growth of the HII regions and the 
last term on the right hand side accounts for "old" reionizations 
due to recombinations inside the HII region. This equation is 
a very rough approximation, as it does not include feedback 
effects, light travel time, and it does not accurately model the 
period when bubbles start overlapping (i.e. Fuu(z) ~ 1). How- 
ever, it can suffice for the crude estimates we are making here. 

In Figure 9, we plot Fuu(z) for several values of M m i n (z) cor- 
responding to redshift-independent values of T vu = 300 K, 10 4 
K, and 10 5 K, from right to left in the figure. The plot assumes 
values of (e„, / esc , N ph/b , C) = (0.1, 0.1, 4000, 4). From the 
figure, it is evident that fairly rapid growth of Fnn(z) is possible 
for Fhii(z) ^0.1. For example, the T vn <; 10 5 K curve goes from 
Fhu ^0.3 to Fhii ~ 0.8 in a the redshift interval z ~ 8.5 — > 7. 
So SNe in roughly 50% of small halos could be wiped out in 
a redshift interval of only Az re <~ 1.5. Comparing with Figures 
4a, 5a or 6, one can see that this could result in a fairly large, 
easily detectable drop in the SNRs. 

To summarize, we outline a likely reionization scenario. An 
early reionization epoch could be driven by minihalos. Feed- 
back processes can stall reionization at a constant or even de- 
creasing filling factor of ionized regions. Then a later popula- 
tion of more massive halos (r v ; r <~ 10 4 - 10 5 K) could complete 
the ionization process on time-scales corresponding to Az le <~ 1 
- 2. The research we present here suggests that this later epoch 
is detectable through an accompanying drop in the SNRs. 
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FIG. 9. — Simple models of the evolution of the filling factor of ion- 
ized regions, Fhii(z), f° r several values of M mm (z) corresponding to redshift- 
independent values of 7/ vlr = 300 K, 10 4 K, and 10 5 K, right to left. 



4. CONCLUSIONS 

We obtain detailed, high-redshift SNe rates which could be 
detected in future infrared surveys, such as JWST. Our method 
uses the extended Press-Schechter formalism to generate high- 
redshift SFRs and SNRs. We empirically calibrate our model 
using the observed properties of local SNe, such as their lightcurves, 
fiducial spectra, as well as dust corrected and uncorrected peak 
magnitude distributions, to predict the detectable high-z SNRs. 

During reionization, the ionizing background radiation heats 
the IGM and could dramatically suppress gas accretion and star 
formation in low-mass halos. As a result, reionization may be 
accompanied by a drop in the SFR and the corresponding SNR. 
The size of the drop is uncertain, since the ability of the gas 
in low-mass halos to cool and self-shield against the ionizing 
radiation is poorly constrained at high redshifts. However, we 
argued, using specific illustrative models, that a relatively sharp 
feature is feasible. Using the SNRs we generate, we then ana- 
lyze the prospects of detecting such a feature with future SNe 
observations, specifically those obtainable with JWST. 

We find that 4-24 SNe may be detectable from z <; 5 at the 
sensitivity of 3 nJy (requiring 10 5 s exposures in a 4.5 band) 
in each ^10 arcmin 2 JWST field. In a hypothetical one year 
survey, we expect to detect up to thousands of SNe per unit red- 
shift at z ~ 6. Our results imply that, for most scenarios, high- 
redshift SNe observations can be used to detect fairly sharp fea- 
tures in the reionization history (occurring over Az re <~ 1 - 3) 
out to z ~ 13, as well as set constraints on the photo-ionization 
heating feedback on low-mass halos at the reionization epoch. 
Specifically, for a wide range of scenarios at z re ^13, the drop 
in the SNR due to reionization can be detected at S/N <; 3 with 
only tens of deep JWST exposures. Our results therefore sug- 
gest that future searches for high-z SNe could be a valuable new 
tool, complementing other techniques, to study the process of 
reionization, as well as the feedback mechanism that regulates 
it. 
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